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^ ■ Abstract: This study explores the reconstruction of a signal using spectral quan- 

tities associated with some self-adjoint realization of an h-dependent Schrodinger 
d 2 

operator — h 2 -^-^ — y(x), h > 0, when the parameter h tends to 0. Theoretical re- 

suits in semi-classical analysis are proved. Some numerical results are also presented. 
ch ! We first consider as a toy model the sech 2 function. Then we study a real signal given 

by arterial blood pressure measurements. This approach seems to be very promising 
in signal analysis. Indeed it provides new spectral quantities that can give relevant 
information on some signals as it is the case for arterial blood pressure signal. 

<N [ 
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cK 1 Introduction 
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o . 

Let y be a positive real valued function on a bounded open interval Q =}a,b[ rep- 
resenting the signal to be analyzed. Following the idea in [9], [10] . in this study we 
interpret the signal y as a multiplication operator, <fi — > y- cj), on some function space. 



5_i ■ The spectrum of a regularized version of this operator, i.e. the Dirichlet realization 



in f2 or the periodic realization of 



H h ,y = -h 2 ^- 2 -y, (1) 
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for a small h, is then used for the analysis of y (see [9], |10]). We will denote by H^ y 
and these realizations and by y one of these two realizations. 
We are interested in analyzing the potential y on a compact K C Vt using the 
negative eigenvalues A^ n , n = 1, • • ■ , and some associated orthonormal basis of 
real eigenfunctions iphm n G N of the Schrodinger operator H£ y . For this purpose, 
we choose A < such that A is a noncritical value of —y, A < inf (-y(a), —y(b))) and 
— y is strictly less than A on K. Then, we show that y(x) can be reconstructed in K 
using the following expression 

Uh(x, A) = -A + 4/i ^ (A - \ hn )^ 4>hn(x) 2 , xeK. (2) 

This expression is different from the one considered in [10] which is given by: 

y h (x,0) =Ah (-Afcn)^ftn(a;) 2 , (3) 

that is corresponding to a problem which is defined, at least at the theoretical level, 
on the whole line, under the assumption that y tends sufficiently rapidly to at oo. 
Indeed, in [10], the whole potential is recovered, using all the negative eigenvalues 
and the associated eigenfunctions of the selfadjoint realization of Hh,y on the line, on 
the basis of a scattering formula due to Deift-Trubowitz |T|. This formula involves 
a remainder whose smallness as h — > is unproved. In addition the numerical 
computations are actually done for a spectral problem in a bounded interval, hence 
it seems more natural to look directly at such a problem. Another advantage of our 
new approach is that we can consider in the same way, with yh(x, A) = y h i(x, A), 

y&, 7 (z,A) = -A+(A (A-A h „r^„(x) 2 )^, xeK. (4) 

7 A hn <A 

for 7 G [0, +oo[ and Uz is a suitable universal semi-classical constant and discuss the 
speed of convergence as h —> of yh^{x, A) in function of 7. 

At the theoretical level, we can either choose the Dirichlet or the periodic realization. 
Agmon's estimates (0 IS]) show indeed that 

y^{x i X)-^r(x,\) = 0(h°°), (5) 

and actually is exponentially small as h — > 0. Nevertheless, at the numerical level, 
the periodic problem seems to give better results in term of accuracy and convergence 
speed. 
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In addition, if we are interested in the analysis of the signal in a specific interval, 
the introduction of a suitable A depending on this interval enables us to have a good 
estimate of this part of the signal with a smaller number of negative eigenvalues. This 
can have interesting applications in signal analysis where sometimes the interest is 
focused on the analysis of a small part of the signal. Let us also mention that 
the interest of our approach is not really in the reconstruction of the signal that 
we already have in fact but in computing some new spectral quantities that provide 
relevant information on the signal. These quantities could be the negative eigenvalues 
or some Riesz means of these eigenvalues. The main application in this study is in 
the arterial blood pressure (ABP) waveform analysis (the signal is then the pressure). 
We refer to [10], |12| where for instance it is shown for A = how these quantities 
permit to discriminate between different pathological or physiological situations and 
also to provide some information on cardiovascular parameters of great interest. 
As described in [10], the parameter h plays an important role in this approach. 
Indeed, as h becomes smaller, the approximation of y(x) by yh(x,0) improves. We 
have the same remark in this study and we will prove the pointwise convergence 
of yh{x,X) (or more generally yh tl {x,X) ) to y(x) when h — > 0. This explains our 
terminology "semi-classical". Note finally that in the applications the choice of h is 
not necessarily very small. Hence the theoretical analysis given in this article has 
only for object to give some information for the choice of an optimal h and possibly 
some A. 



Let us now present our main result which will be later obtained as a particular case 
of a more general theorem. 

Theorem 2.1. 

Let y be a real valued C°° function on a bounded open set Q. Then, for any pair 
(K, A) such that K is compact and 



2 Main results 





and, uniformly for x G K , 



we have 



y(x) 



A + lim4/i V (A - \ hn )z ^ hn {xf 




(7) 
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where Xh n and iph n denote the eigenvalues and some associated L 2 -normalized real 
eigenfunctions of H£ y in Q. Here we recall that H£ y denotes either H® y or H^ y . 

Remark 2.2. 

We do not impose in the statement of Theorem \2.1\ the positivity of y and the neg- 
ativity of X. If we want to recognize the signal y on a given compact K , then any 
choice of A such that ([6]) is satisfied is possible. There are at this stage two con- 
trary observations. To choose A as small as possible will give the advantage that we 
need less eigenvalues and eigenfunctions to compute, but this is only true in the limit 
h — y and this imposes to compute more eigenvalues! For a given h, take a larger A 
seems numerically better. An explanation could come from another asymptotic anal- 
ysis which will not be done here. Hence our theorem gives just some light on what 
should be done numerically and one has to be careful with their interpretation. 

We will obtain the proof by using a suitable extension of Karadzhov's theorem on 
the spectral function. One can suspect that the convergence could be better due to 
the fact that A is not critical, but it should not affect so much the computation if 
y(K) + A is sufficiently large. 
So the basic idea is that 

h J2 X(A - X hn ) ^hn(x) 2 = (27T)- 1 I X (X + y(x) - e) <% + 0(h) , 

uniformly with respect to x in a compact and that we can have a complete expansion 
in h if supp \ CC]0, +oo[. In what follows we omit the subscript y in Hh y . The main 
point here (in the case of the line) is, following [3], that x(A — Hh) can be considered 
as an /i-pseudodifferential operator whose /i-symbol admits an expansion 

q(x,{,h)~J2h j Qj(x,0 (8) 

where the qj{x,£) have compact support, 

qo(x,0 = x(X + y(x)-e) (9) 

and 

qi(x,O = 0. (10) 
This means that, for u G C^°(IR), we can write 

( X (A - H h ) u) (x) = (27T/*)- 1 J ^ (»-»)-€ ff (£±^, e, h)u{ y )dydi . (11) 
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We are then considering the restriction to the diagonal of the distribution kernel Kh 
of X (A - H h ) 

K h {x,y) := (27rh)' 1 J e iM<qfi±V,Z,h)d£, (12) 
which at (x, x) becomes 

K h (x, x) := (21th)- 1 j q{x, e, /iR , (13) 
and admits the expansion in powers of h 

K h (x,x) ~ (2Trh)- 1 Y / h j j qsMdt. (14) 

Moreover, when integrating over x, we get the trace of %(A — H^) 

Tr x(A - fr h ) = (27T/1)" 1 y g(x, £, /i)da;d£ ■ (15) 
We also note for future use (see (2.3) in j3]) 

q 2 (x, = ^x"(A + y(z) - £ V(*) + ^X'"(A + y{x) - fXV'Wf ~ y'{xf) . (16) 

The hope is that when x is replaced by A+, with 7 > 0, we keep a remainder in 
0(h l+1 ) in the right hand side. This reads 

h £(A - A, n )X Vhn(x) 2 = (2-kY 1 j{\-e + y(x))l<% + 0(h 1+ ^) . (17) 

1 

This suggests the change of variable £ = r/(A + y(x))+ and leads to 

h ^(A ~ X hn)l *hn(x) 2 = (A + y(x))r i (27r)- 1 c 7 + 0(h 1+ ~<) . (18) 

with 

c 7 := J (1 - V 2 )ld V . (19) 

The hope is that the remainder will be uniform for x in a compact X such that 
y(K) C ]-A,+oo[. 

This theorem is proven in [S] (with some details given in [7J) for 7 = 0, see in the 

next section. For our application, we need to have at least the case 7 = - but we 
will also suggest in our conclusion that considering 7 larger could be better. 
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3 Former results 

We recall one of the basic results of Karadzhov [7|. We come back to the more 
standard notation by writing V = —y. 

Theorem 3.1. 

d 2 

Let Ah = —h 2 —— + V, with a C°° potential V on the line with V tending to +00. Let 
dx l 

eh be the spectral junction associated to Ah such that e^(A, x, y) = iphj{x) V'h •(?/)■ 

Xj<X 

Then, for any compact K such that K C V r ~ 1 (] — 00, A[), we have, for x e K 

e h (\,x,x) = n~ 1 (\-V(x))h~ 1 + 0(l), h -> 0, (20) 

uniformly in K . 

The next theorem was established by Helffer-Robert jl] in connection with the anal- 
ysis of the Lieb-Thirring conjecture. This will not be enough because we will need 
a point-wise estimate and this is an integrated version, but this indicates in which 
direction we want to go. 

Theorem 3.2. 

Let V G C°°(R), with 

—00 < inf V < lim inf V . 

\x\ — >+oo 

Let A e] inf y, liminf^i^+oo y[ and suppose that A is not a critical value for V. We 
denote by : 

S,(h,X)= (A-A, n ) 7 , 7 >0 (21) 

the Riesz means of the eigenvalues Xh n less than A of Ah- Then for < 7 < 1, we 
have: 

S 7 (h, A) = - \L* J (A - V(x))l + * dx + O(h^) (22) 
where | • | + is the positive part and Uz, known as the classical Weyl constant, is given 

by 

(23) 



. d r( 7 + i) 

^ 20Fr(7-^ 



2 - 



3 1 
Note that = c 7 (27r)~ 1 and we recall that r(— ) = . Hence in particular Lq = — 

and Lf 
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4 Proof of the main theorem 



4.1 Pointwise asymptotics for the Riesz means 

The statement of the main theorem will correspond to the case 7 = — of the following 
Theorem 4.1. 

d 2 

Let Hi be the realization of Hh = —h 2 -—r + V, with a C°° potential V = —y (We 

ax 2 

can either consider or H^ er ). Let el be defined by 

el{X,x,y) = ( A - X hjyiphj(x)ip hj (y) . 
x hj <x 

For any pair (K, A) satisfying (J6]) , then we have, for x e K, 

el(X,x,x) = (2tt)- 1 (A- V{x)Y f+ ^o 1 hr 1 + 0(W), h^O, (24) 
uniformly in K . 



We refer to the heuristics starting from (117j) for understanding the main term in (j2] 
As in the proof by Helffer- Robert of Theorem 13.21 [4], we can distinguish two steps 
corresponding to the contribution which is close to A and to the contribution which 
is "far" from A. This will be done by a cut-off in energy. 

4.2 Far from A 

We consider a function \ with suppx CC] — 00, 0[ and consider the expression 

fl x (X,x,x) := J^xiK - A)|A - X hj \^ hj {xf , (25) 
3 

and we prove that 

Proposition 4.2. 

If (A, K) satisfies the condition 

X<mi(V(a),V(b)), (26) 
and if K is compact in }a,b[ then, uniformly on K, we have 



ft A\,x,x) 



(a*)" 1 /!- 1 [Y, a U x )A ■ ( 27 ) 
Wo / 
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where 

°4x(s) = / x(A-e 2 -^))iA-nx)-e 2 r^. (28) 

and the aj x are C°° functions on ]a, b[. 

4.3 Close to A 

We now consider 

gl x (X,x,x) := (1 - X(A« - A))|A - X hj \^ hj {xf • (29) 

This quantity involves only the eigenvalues (and corresponding eigenfunctions) which 
are close to A. Of course, we have 

el(X, x, x) = fl x (X, x, x) + gl tX (X, x, x) , (30) 

and for coming back to our main theorem we keep in mind that 

e° h (X,x,x) = e h (X,x,x) . (31) 

Assumption (126|) together with Agmon estimates permits (in the two cases, Dirichlet 
or periodic) to reduce to a global problem on R where V is replaced by a new 
potential V coinciding with V on K. Moreover, we can impose conditions on V 
at oo permitting to use the global calculus of Helffer-Robert calculus (or the semi- 
classical Weyl calculus) and to get the result. 

The support of x is now chosen so that % = 1 on ] — oo, — e] with e > small enough. 
The second proposition is 

Proposition 4.3. 

For any pair (K, X) satisfying Q and for a sufficiently small e > 0, we have uni- 
formly for x G K 

gl a (X, x, x) = (2ny 1 h- 1 {fijx) + 0(h)) , (32) 

where 

nq x (x) = J(i-x(x-e-v(x)))(x-v(x)-e)idn. m 
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Here the proof is a rather direct consequence of Karadzhov's paper. We are indeed 
analyzing 

gl >x (X, x,x)= (1 - x(A - n)) (A - /i) 7 de h (fi, x, x) . 
J fi<\ 

Using an integration by parts, we obtain 



g% (\,x,x)= / ^(A- fi)e h (fi,x,x)dfi, 
J n<x 

with 

m = (a - xmir . 

We can then apply the estimate (1241) we had for eh(fi,x,x) which is uniform for \i 
close to A and we obtain by integration over ] — oo, A] (and the reverse integration 
by parts) the result modulo 0(h). 

Theorem 4.4. 

For any pair (K, A) satisfying ([6]) and sufficiently small e > 0, we have uniformly 
for x G K 



(2n)- 1 h- 1 ( Pljx)h' + 0(h^)Y (34) 

\0<£<7+l / 



with 



and 



/3 y*) = j (i- X (\-e-v(x)))(\-v(x)-ey + dn, (35) 

el(X,x,x) = (2n)- 1 h- 1 ( ^ aj(x)h e + 0(h 1+ ^) J . (36) 

\0<£<7+l / 

with 

al{x) = f(\-V(x)-ey + dt. (37) 
Here the j3j x and a] are C°° functions in a neighborhood of K. 

A direct proof should be given for having this remainder. One should follow Karadzhov's 
proof (in its easy part because we are far from V(x) = A) and improve the Tauberian 
theorem used in this paper, using the ideas of jl] (in a non integrated version). It 
is also based on the approximation of (1 — x)Hh exp(— it^) by a Fourier Integral 
Operator. We do not give the details here. 
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More information on the coefficients 

We get from p[| that ol[{x) = . Actually only the integrated version of this claim 
is given (see (0.12) there and have in mind that the subprincipal symbol vanishes) 
but coming back to the proof of [3] gives the statement . It is also proven by Helffer- 
Robert |3] that is not identically zero. The computation is easier for 7 > 2. We 
follow what was done in [4] (Equation 2.20), but we can no more integrate in the x 
variable, so the simplification obtained in this paper by performing an integration 
by parts in the x variable is not possible. We get, for 7 > 2, (see (fIE]l ) 

a$(x) = L 7 (A - V(x)y~^V"(x) + L'^X - V(x)y~^V \x)\ (38) 

where L 7 and L 7 are universal computable constants. So when increasing 7, we 
suspect that we improve the semi-classical approach of y(x) when we increase 7 from 
to 1 and that we do not improve anymore for larger 7. Fig. [7] confirms this guess 
for 7 = 2. The asymptotic behavior of y(x) — yh n {x) is indeed given by h 2 5l(x), 
where 5 2 is a C°° function in a neighborhood of K directly computable from a 2 {x) 

and V. In the case 7 = -, we see on the contrary some oscillation. 



5 Some numerical examples 

To study the "concrete" validity of our formula, we have performed some numerical 
tests using MATLAB software. We have chosen to use a pseudo-spectral Fourier 
method instead of a finite differences method for the discretization of the prob- 
lem. Indeed, Fourier method gives better results in term of accuracy and speed 
convergence. However Fourier method requires periodic boundary conditions so the 
numerical tests have been done on the periodic realization H^. 
We consider a grid of M equidistant points xj, j = 1, • • • , M such that 

a = X\ < x 2 < • • • < xm-i < %m = b. (39) 

We denote yj and i/jj the values of y and at the grid points Xj, j — 1, • • • , M 

y 3 = y{ Xj ), ^ = ii{xj), j = 1, ■ • • , M. (40) 

Therefore, the discretization of leads to the following eigenvalue matrix problem 

(-h 2 D 2 - diag (Y)) ± = -X±, (41) 

where diag (V) is a diagonal matrix whose elements are yj, j — 1, • • • , M and if) = 

["01 t l J 2- l "•• tpM-i 4'm} T - D 2 is the second order differentiation matrix for a pseudo- 
spectral Fourier method [TU], [15] . Note that the analysis is only relevant if M is 
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large, and h can not be too small in comparison with the distance — between 

two consecutive points. We denote Nh the number of negative eigenvalues of H^y 
and Nh t \ the number of negative eigenvalues less than A. 

Different values of the parameter h, A and 7 have been considered, sometimes outside 
the probable domain of validity of the theoretical analysis. We concentrate our 
analysis on two examples. We first look at the sech 2 function which is a regular 
function and then we consider the case of arterial blood pressure signal. The latter 
is given by measured data at the finger level provided by physicians and it is difficult 
in this case to speak of regularity of the signal. 

5.1 The sech 2 signal. 

We first study as toy model the example of the sech 2 function defined on [0, 10] by 

y(x) = sech 2 (x-5). (42) 

It is a well studied potential when considered on the whole line with explicitly known 
negative spectrum for the associated Schrodinger operator for some specific values 
of h. 

We analyze the reconstruction of a part of the signal given by — y < Ai, with Ai = 
—0.8. For this purpose, we use Formula ([3]) and compare the reconstruction of y(x) 
with this formula with its reconstruction with yh n {x, A) for different values of h and 

A. Fig. [1]- fig. H] illustrate the results for 7 = -. We notice that for h = 0.01, 

the reconstruction of this part of the signal is better with y^(x,0). However, when 
h = 0.001 the reconstruction is better with yh(x, A). 

Then, we analyze the error with different values of 7 and h, A being fixed. Fig. [5]- 
fig. [7] illustrate the results for A = —0.5. The optimal choice of 7 for fixed h and A 
seems to be 7 = 1. Note also that the error for 7 = 2 is regular as it was explained 
previously. 

5.2 The arterial blood pressure (ABP) signal 

We analyze in this section the results obtained when the signal y(x) is an ABP 
signal. ABP plays an important role in the cardiovascular system and is used in 
clinical practice for monitoring purposes. However the interpretation of ABP signals 
is still restricted to the interpretation of the maximal and the minimal values called 
respectively the systolic and diastolic pressures. No information on the instantaneous 
variability of the pressure is considered. Recent studies have proposed to exploit the 
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h=0.01, *.=- 0.1, N =99, N. 




h=O.0O1, ^.=-0.1, N. =991, N,,,=684 



(a) h = 0.01 



0.024 
0.022 
0.02 
0.018 
i 0.016 




4.5 4.6 4.7 4.8 4.9 5 5.1 5.2 5.3 5.4 5.5 



(b) h = 0.001 



Fig. 1: Relative error with A = and A = —0.1 for the sech 2 example 



h=0.01, X=-0.3, N. =99, N. 



_ _ _ Error with y (x,X) 



h=0.001, X--0.3, N =991, N =452 





4.5 4.6 4.7 



4.9 5 5.1 5.2 5.3 5.4 5.5 



(a) h = 0.01 (b) h = 0.001 

Fig. 2: Relative error for A = and A = —0.3 for the sech 2 example 
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h=0.01, 3w=-0-5, N. =99, N hx =29 
_ _ m Error with y.(x,X) 



h=0.001, X--0.5, N, =991, K , -293 





4.5 4.6 4.7 4.8 4.9 5 5.1 5.2 5.3 5.4 5.5 



(a) h = 0.01 (b) h = 0.001 

Fig. 3: Relative error for A = and A = —0.5 for the sech 2 example 



h=0.01, ?.=-0.7, N,, =99, r>1 . 



h=0.001, X=-0.7, N =991, N =163 
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_ _ _ Error with y (x,l) 
Error with y (x,0) 



_ Error with y (x.,X) 
_ Error with y.(x,0) 




4.5 4.6 4.7 4.8 4.9 5 5.1 5.2 5.3 5.4 5.5 



(a) h = 0.01 (b) h = 0.001 

Fig. 4: Relative error for A = and A = —0.7 for the sech 2 example 
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h=0.1, i.=-0.5, N, =10, N. . =3 



h=0.1, ?.=-0.5, N h =10, N, . =3 
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4.5 4.6 4.7 4.8 4.9 5 5.1 5.2 5.3 5.4 5.5 



(a) 



(b) 



Fig. 5: Reconstruction of a part of the sech 2 signal (a) and relative error (b) for 
different values of 7 and h — 0.1 



h=0.01, a.=-0.5, N. =99, N,,,=29 




h=0.01, X=-0.5, N. =99, N,,,=29 




4.5 4.6 4.7 4.8 4.9 5 5.1 5.2 5.3 5.4 5.5 



(a) 



(b) 



Fig. 6: Reconstruction of a part of the sech 2 signal (a) and relative error (b) for 
different values of 7 and h = 0.01 
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4.4 4.6 4.8 5 5.2 5.4 5.6 4.5 4.6 4.7 4.8 4.9 5 5.1 5.2 5.3 5.4 5.5 



(a) (b) 

Fig. 7: Reconstruction of a part of the sech 2 signal (a) and relative error (b) for 
different values of 7 and h = 0.001 

ABP waveform in clinical practice by analyzing the signal with a semi-classical signal 
analysis approach (see for example |10| . |11|). The latter proposes to reconstruct the 
signal with formula ([3]), corresponding to the case A = in order to extract some 
spectral quantities that provide an interesting information on that signal. These 
quantities are the eigenvalues and some Riesz means of these eigenvalues. These 
quantities enable for example the discrimination between different pathological and 
physiological situations [12] and also provide information on some cardiovascular 
parameters of great interest as for example the stroke volume [11]. We are interested 
in this study in the reconstruction of a small part of an ABP pressure beat illustrated 
in fig. [8] with formula @. We consider different values of h, A and 7 as it is described 
in fig. [9] - fig. [22] where — y is represented. In fig. [9] - fig. [16] a zoom on the 
reconstructed part of the signal is represented. Note that the signal in this case 
is not a regular function but only known from data measurements for a sequence 

Xj — a + (j — 1)tt — -r U = !)•••) M) for some integer M. We recall that in our 

application the x variable represents the time. The time between two consecutive 
measurements is 10 -3 seconds. 

Fig. |9] and fig. [10] as well as fig. [T7] and fig. [18] clearly show that, as h decreases, 
the approximation improves. Fig. [TT] and fig. Q3] suggest that, for h fixed, as the 
distance of A to —y(K) increases, the estimate improves in K. Fig. [T2] fig. EH fig. 
dU fig. [inland fig. [2U] illustrate the influence of 7. The differences between the errors 
of reconstruction in the three cases are not significant. However, the error for 7 = 2 
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is not regular as it was the case for the toy model. This can be explained by the 
non regularity of our signal. Note also that these reconstruction formulae appear as 
a regularization of the signal. 

Fig. [151 fig. [161 fig- ED and fig. [22] compare the case A = to the case A = — 70 for 
the first example and to the case A = —65 for the second example. The errors are of 
the same order of magnitude. So it confirms our idea that we have not to take all the 
eigenvalues to estimate the signal. Indeed, we have found in the two smaller 
A such that the estimate in K is satisfactory with a smaller number of eigenvalues 
and corresponding eigenfunctions to use. 

Fig. [23] is a good illustration of many semi-classical properties. First we see that, 
below the energy —60, there is a concentration of the eigenvalues near the critical 
points of —y. Far above —60, the behavior is asymptotically given (for fixed h) by 

d 2 

the eigenvalues of the periodic realization of — h 2 — - in la, b[ with a and b describing 

ax 2 

the beginning and the end of an arterial blood pressure beat respectively. 

6 Conclusion 

In continuation of [5J, |10| . |11| . we have explored the possibility of reconstructing 
the signal y using spectral quantities associated with some self-adjoint realization 
of an /i-dependent Schrodinger operator (h > 0) —h 2 -^ — y(x), the parameter h 
tending to 0. Using on one hand theoretical results in semi-classical analysis and on 
the other hand numerical computations, we can formulate the following remarks. 

• As h (h > 0) decreases (semi-classical regime), the approximation of y by 
y/ lj7 (x, A) improves. 

• Semi-classical analysis suggests also to take 7 > 1. 

• However, the number of negative eigenvalues to be computed increases as h 
becomes smaller. So the numerical computations become difficult if h is too 
small. 

• For a given interval K, a clever choice of A makes possible to get a good 
approximation of the signal in K in the semi-classical regime with a smaller 
number of eigenvalues. 

• Finally, the numerical computations suggest also (but we are outside the the- 
oretical considerations of our paper) that one can consider a choice of h (not 
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necessarily small) and hope a good reconstruction of the signal by choosing 
appropriate values of A and 7. This should be the object of another work. 
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(a) Example 1 (b) Example 2 

Fig. 8: An example of an arterial blood pressure beat measured at the finger. The 
part of the signal that we want to reconstruct is in red. 
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h=0.1, N =31, N hx =9, ¥*1/2 



h=0.1, N =31, N =9, y=1/2 




- Error with k=-70 

- Error with X=0 




0.09 0.1 0.11 0.12 0.13 0.14 0.15 



(a) 



(b) 



Fig. 15: Reconstruction of a part of the ABP signal (a) and relative error (b) for 
A = and -70, h = 0.1 and 7 = \ 



h=0.01, N =305, N =86, y=1/2 



h=0.01, N =305, N =86, y=1/2 




- Error with X=-70 

- Error with ^=0 
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(a) 



(b) 



Fig. 16: Reconstruction of a part of the ABP signal (a) and relative error (b) for 
A = and -70, h = 0.01 and 7 = \ 
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Fig. 19: Reconstruction of a part of the ABP signal (a) and relative error (b) for 
A = -65, h = 0.1 and 7 = §,1,2 
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h=0.1, N.=31, X=-70 

n 

Or 
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-140 1 1 1 1 1 1 1 ' 

1.6 1.8 2 2.2 2.4 2.6 2.8 3 

t 

Fig. 23: The negative eigenvalues of the Schrodinger operator with a potential given 
by ABP signal 
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